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Abstract 


A description  is  given  of  the  design,  operation  and  performance  of 
a radically  new  type  of  sonic-boom  generating  flap- valve.  Included  are  the 
ancillary  cam,  clutch,  flywheel  and  electric  motor  system,  as  well  as  a much 
larger  air  reservoir  than  used  previously.  An  rqpdated  and  greatly  extended 
analysis  describes  the  time  varying  reservoir  conditions  (e.g.,  pressure)  and 
mass-flow  rate  of  air  throvigh  the  flap  valve,  as  well  as  the  wave  motion  or 
characteristics  of  the  travelling  N-wave  in  the  pyramidal  concrete  horn  (25  m 
long  with  a 3 m X 3 ® base)  of  the  UTIAS  travelling-wave  sonic-boom  simulation 
facility. 
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1.  IMIRCDUCTION 


With  each  new  year  more  of  mankind  1b  being  sxibjected  to  additional 
noise  and  its  consequent  increased  stress.  Research  into  the  effects  of  noise 
on  humans,  animals  and  structures,  as  well  as  its  abatement,  is  also  expanding 
at  an  accelerated  pace,  fortunately.  Results  of  such  work  can  provide  a base 
for  establishing  safe  guidelines  of  noise  exposure.  One  inportant  and  currently 
active  area  of  noise  research  is  the  study  of  the  effects  of  inpiilse  sound. 
Furthermore,  a particularly  inportant  inpiolse  sound  is  the  sonic  boom,  whose 
inpact  is  now  being  assessed  before  supersonic  transport  (SST)  aircraft  such  as 
the  Anglo-French  Coii-orde  and  the  Soviet  TU-144  are  introduced  into  extensive 
commercial  service. 

A sketch  of  the  shock-wave  pattern  emanating  from  a SST  aircraft  and 
the  trailing  N- shaped  overpressure  signatiure  at  the  Earth's  surface  are  shown 
in  Fig.  1 (Ref.  1).  Furthermore,  a better  illustration  of  the  groiuid-intercepted 
sonic-boom  path  appears  in  Fig.  2.  The  most  inportant  parameters  which  are 
commonly  used  to  describe  the  various  parts  of  the  signature  include  the  peak 
overpressiire,  rise- time,  duration  (or  wave  length) , and  wave  form  which  may 
deviate  somewhat  from  the  ideal  N-shape.  It  is  worth  notix:g  that  respective 
values  of  peak  overpressure,  duration  and  rise- time  are  100  N/i^,  300  ms  and 
1 ms  for  a typical  sonic  boom  from  a ciurrent  SST  and  also  from  a large  military 
sipersonic  bonber.  In  the  case  of  a shorter  supersonic  fighter,  only  the 
duration  is  significantly  different,  being  correspondingly  shorter  at  about 
100  ms. 


In  order  to  successfully  simulate  a sonic  boom  from  a SST  aircraft, 
any  simulation  facility  must  have  the  capability  of  repeatedly  producing  an 
N-wave  having  the  correct  peak  oveipressure,  rise- time  and  duration.  In  the 
late  1960's  and  early  1970' s,  many  different  types  of  sonic-boom  simulation 
devices  and  techniques  were  proposed  and  nu)st  simulators  or  prototypes  were 
constructed  and  tested.  A description  of  such  simulators  and  theii  test  results 
cannot  be  given  here,  but  the  interested  reader  can  find  much  information  in 
various  review  articles  or  reports  (Refs.  2,  3 and  4).  However,  it  is  worth 
mentioning  that  most  simulators  were  only  partially  successful  (e.g.,  produced 
an  N-wave  having  the  correct  peak  overpressure  and  rise  time  but  not  duration), 
and  only  a few  met  all  of  the  desired  requirements  adequately  (Refs.  5 to  10) . 

The  sonic-boom  program  at  UTIAS  started  in  the  late  1960's.  For 
exeaple,  see  Ref.  1 for  a review  of  the  work  conpleted  up  to  1974.  A good  part 
of  this  program  was  the  development  of  suitable  sonic-boom  simulators.  One 
portable  sonic-boom  simulator  in  the  form  of  a shock  tube  (11  kg,  1 m long) , 
having  a constant-axea  driver  and  exponential  horn,  can  easily  be  transported 
and  operated  by  one  person  to  conduct  wildlife  field  tests  (Refs.  1 and  11). 

A simulated  short-djiration  sonic  boom  can  be  produced  and  directed  at  wildlife 
in  its  natural  habitat,  in  order  to  study  its  startle  response  and  subsequent 
behaviour. 


The  first  of  two  major  laboratory  facilities  at  UTIAS  is  the  Loud- 
speaker-Driven Booth,  which  can  easily  accommodate  one  hvnnan  subject  or  small 
caged  animals  in  its  solidly  built  and  sealed  chamber  (about  2 n^)  to  facilitate 
response  tests  to  a simulated  full-scale  sonic  boom  (Refs.  1,  9 and  12).  The 
second  major  laboratory  facility  is  the  Travelling-Wave  Horn  (Refs.  1,  4,  10,  12, 
13  and  l4) , and  its  essential  features  are  shown  in  the  elevation  and  plan  views 
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of  Fig.  3-  The  large  horizontal  pyramid  is  25  m long,  has  a 3 m square  base, 
and  is  made  mainly  of  concrete.  Near  the  horn  apex  a sonic-boom  generator  in 
the  form  of  either  a mass-flow  valve  or  a shock-t\jbe  driver  is  used  to  control 
the  air  flow  from  the  high-pressure  reservoir  into  the  horn.  The  expanding  air 
flow  near  the  horn  apex  simulates  the  expansion  process  occurrir.g  in  a weak 
spherical  explosion  and  produces  the  travelling  N-wave  or  simulated  sonic  boom 
which  propagates  from  the  small  to  the  large  end  of  the  horn.  For  studies  of 
the  effects  of  sonic  boom  on  humans  and  animals,  a human  subject  or  small  caged 
animals  can  be  put  in  the  interior  test  section  or  alternatively  in  the  psycho- 
acoustic test  room  wlx.  _h  is  joined  to  the  horn  interior  by  a suitably  sized  open 
window.  For  structural  tests,  one  large  or  many  small  panels  representative  of 
a house  interior  wall  can  be  installed  in  the  cutout,  or  one  or  more  walls  of 
the  psychoacoustic  best  room  can  be  used  directly  for  full-scale  wall  response 
and  fatigue  tests.  Many  such  studies  and  also  others  which  utilized  the 
Travelling-Wave  Horn  have  already  been  conpleted  (see  Refs.  13  and  15  to  22). 

In  the  shock-tube  mode  of  operation,  high-pressure  air  contained  in  a 
constant-area  or  pyramidal  driver  is  suddenly  released  by  breaking  a diaphragm 
which  initially  separates  the  driver  gas  from  the  lower  pressure  air  in  the 
horn.  The  expanding  driver  gas  generates  the  travelling  N-wave  or  simulated 
sonic  boom.  Althoxigh  the  peak  overpressure  of  the  N-wave  can  be  easily  controlled 
to  be  less  than,  equivalent  to  or  greater  than  that  of  an  actual  sonic  boom,  the 
duration  is  short  (x;p  to  20  ms)  and  the  rise  time  is  rapid  ( about  20  jis) . Such  a 
short  N-wave  can  be  useful  for  certain  human  and  animal  response  studies,  but  it 
is  more  suitable  for  studies  of  N-wave  diffraction  over  and  into  a model  building, 
and  also  N-wave  propagation  over  a reduced-scale  topology. 

For  the  mass- flow- valve  mode  of  operation  of  the  Travelling-Wave  Horn 
a specially  designed,  pneumatically-operated  plug  valve  is  used  to  control  the 
mass-flow  rate  of  air  from  the  high-pressure  reservoir  into  the  horn  (Figs.  4 
and  5)  • The  peak  oveipressure  and  wave  length  of  the  simxalated  sonic  boom  can 
be  controlled  independently  such  that  either  one  is  less  than,  equivalent  to  or 
greater  than  that  of  an  actual  sonic  boom.  The  rise  time  of  the  simvilated  sonic 
boom  is  typically  3 to  6 ms.  The  interior  of  the  horn  is  equipped  with  a special 
high-frequency  sound  absorber  or  low-pass  acoustic  filter  for  removing  jet  noise 
from  the  passing  N-wave  (Figs.  3 and  6).  This  undesirable  jet  noise  is  generated 
by  the  high-speed  turbulent  air  flow  at  the  plug  valve,  and  it  is  superposed  on 
the  simulated  sonic  boom.  The  open  base  of  the  pyramidal  horn  is  covered  by  a 
specially-designed  reflection-eliminator  (Fig.  3)j  which  is  in  the  form  of  a 
recoiling  porous  piston.  This  device  adequately  minimizes  the  undesirable  re- 
flected wave  which  is  generated  when  the  simulated  sonic  boom  reaches  the  large 
end  of  the  horn. 

Jet  noise  superposed  on  simulated  sonic  booms  generated  with  the  plug 
valve  is  well  illustrated  in  Figs.  7 and  8 for  different  N-wave  ainplitudes  and 
durations  (Ref.  10)  Without  the  use  of  the  jet-nbise  absorber  (first  column) 
the  superposed  noise  is  normally  very  severe,  whereas  with  the  absorber  (second 
column)  the  noise  is  markedly  less  intense.  However,  even  in  the  latter  case, 
the  subjective  loudness  of  the  boom  can  be  increased  siibstantially  by  the  presence 
of  the  jet  noise,  and  thus  results  of  human  and  animal- response  tests  can  be 
affected  (Ref,  13).  Fvirthermore , since  structural  panels  respond  not  only  to  the 
basic  N- shaped  pressure  signature  but  also  to  the  lower  frequencies  of  the  jet 
noise,  panel  response  is  affected  by  the  jet  noise  (Refs.  l8  and  19)*  Consequently, 
it  was  concltided  in  Ref.  10  that  the  Travelling- Wave  Horn,  when  operated  with  the 
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plug  valve  and  Jet-noise  absorber,  is  limited  in  performance  by  the  Jet  noise 
to  N-waves  having  a peak  a2i?)litude  less  than  2CX)  N/m?  and  a duration  shorter 
them  1^0  ms. 

The  possibility  of  greatly  inproving  the  performance  of  the  facility 
was  pointed  out  in  Ref.  10.  An  analysis  of  the  Jet  noise  showed  that  the  noise 
could  be  reduced  to  inconsequential  significance  by  using  a valve  which  has  a 
much  larger  throat  area  (e.g.,  tenfold  larger  at  about  350  cm2).  For  the  same 
mans-flow  rate  of  air  through  the  valve,  which  would  produce  an  equivalent 
anplitude  and  dvtration  N-wave,  the  larger  valve  would  result  in  a much  lower 
flow  speed  at  the  valve  throat  (well  below  sonic  speed)  and  thus  markedly- 
less  intense  Jet  noise.  It  was  realized  that  scaling  up  the  plug  valve  to  a 
much  larger  size  might  be  inpractical.  The  resulting  larger  and  more  massive 
plug  might  be  very  difficult  to  accelerate  quickly  to  high  speed  when  needed, 
maintain  a constant  high  speed  and  decelerate  in  order  to  control  with  precision 
the  mass-flow  rate  of  air  through  the  valve.  Consequently,  a radically  new  mass- 
flow  valve  was  designed  and  constructed,  and  the  Travelling-Wave  Horn  was  modified 
to  operate  with  this  large  flap  valve. 

A detailed  description  is  given  of  the  design,  operation  and  performance 
of  the  UTIAS  Travelling-Wave  Horn  with  the  new  flap  valve.  This  description 
covers  the  basic  elements  of  the  facility  including  the  air-conpressor  and  large- 
reservoir  system,  pyramidal  horn,  reflection  eliminator,  and  flap  valve  with  its 
unique  electric-motor  drive,  flywheel,  fast-acting  clutch  and  cam  system.  This 
is  followed  by  an  ipdated  and  greatly  extended  analysis  to  describe  the  time- 
varying  reservoir  conditions  (e.g.,  reservoir  pressure),  valve  operation  and 
wave  motion  in  the  pyramidal  horn.  The  arplitude,  duration  and  wave  form  of 
the  simulated  sonic  boom  can  be  predicted  successfully. 


V.  2.  DESCRIPTION  OF  THE  TRAVELLIHG-WAVE  HORN 

'i  j 

2.1  Pyramidal  Horn 

cj  A plan  view  of  the  UEIAS  sonic-boom  laboratory  and  Travelling-Wave 

r j Horn  when  operated  with  the  flap  valve  is  shown  in  Fig.  9.  The  pyramidal  horn 

is  enclosed  at  the  small  end  by  a building  called  the  control  room  and  at  the  ; 

1 large  end  by  the  test  room,  as  depicted  in  Fig.  9 and  also  shown  in  Fig.  10. 

Pictvtres  of  those  parts  of  the  horn  contained  in  the  control  suid  test  rooms 
appear  in  Figs.  11  and  12,  respectively.  The  interior  of  the  horn,  looking  from 
j the  base  towards  the  apex,  is  shown  in  Fig.  13,  and  the  opposite  view  of  the 

interior  appears  in  Fig,  l4.  The  illusory  effects  that  the  horn  appesirs  infi- 
nitely long  in  Fig.  13  and  has  no  divergence  in  Fig.  l4  are  also  experienced  \ 

when  one  stands  inside  the  horn.  j 

The  pyramided  horn  is  25  m long,  has  a square  base  which  is  3 ® on  each  / j 

side,  and  has  a total  divergence  angle  of  7.2  degrees.  The  first  part  of  the  j 

j horn  which  is  2.4  m long  is  made  of  2.5-cm-thick  steel  plate,  and  it  is  supported 

by  a special  stand  (Fig.  11).  The  remainder  of  the  horn  is  made  of  steel- 
I reinforced  concrete,  and  this  monolithic  structure  also  has  extremely  rigid  walls 

which  are  20  cm  thick.  The  steel  and  concrete  portions  of  the  horn  were  purposely 
j designed  to  have  very  rigid,  nonporous  and  flat  walls,  in  order  to  minimize 

undesirable  wave-energy  losses  and  resulting  wave- form  distortion  as  the  simulated  i 

sonic  boom  propagates  from  the  small  to  the  large  end  of  the  horn.  The  horn  walls 
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are  plane  to  within  0.2^  of  the  duct  width  where  the  wall  perturbation  occurs, 
and  such  protuberances  shotild  not  produce  significant  transverse  waves  in  the 
horn. 

The  cross-sectional  area  of  the  horn  (A)  increases  continuously  with 
radial  distance  (r)  measured  from  the  projected  horn  apex,  as  given  by  the 
following  expression; 


A = Ag(r/r^)^  = (r/8)^ 


(1) 


f 
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The  area  at  the  large  end  (Ae)  equals  9*30  and  the  corresponding  radius 

(^e)  equals  24.38  m.  This  expression  for  the  plane  cross-sectional  area  is  a 
good  approximation  for  the  curved  surface  area  associated  with  the  spherical 
wave  front  of  the  simulated  sonic  boom,  and  it  is  used  for  sinplicity  in  the 
analytical  work.  Note  that  the  percentage  difference  between  the  plane  and 
curved  surface  areas  amounts  to  only  0.3^  for  the  OTIAS  horn,  because  the 
divergence  angle  of  7.2  degrees  is  relatively  small. 

A pyramidal  horn  instead  of  a hyperbolic  or  wedgy  duct  was  selected 
because  this  type  of  horn  forms  a solid  angle  of  a sphere.  Consequently,  the 
expansion  flow  process  of  a weak  spherical  explosion  (Ref.  l4)  can  be  simulated 
and  a travelling  N-wave  produced  from  the  onset  near  the  horn  apex.  Only  a small 
portion  of  a full  sphere  is  incorporated  in  the  horn  in  order  to  minimize  the 
source  energy  required  to  produce  the  simulated  sonic  boom.  The  source  energy 
required  for  the  operation  of  the  ITEIAS  facility  is  about  three  orders  of 
magnitude  (l/2^6ir)  less  than  that  needed  for  the  case  of  a full  sphere,  >Ailch 
is  calciilated  by  taking  the  ratio  of  the  cross-sectional  area  of  the  horn 
(r2/64)  to  the  fxd.1  area  of  the  corresponding  sphere  (477r2) . The  source  energy 
cannot  be  reduced  much  further  because  a sufficiently  large  interior  test  section 
is  required  (about  2.5  m square)  to  accommodate  a hianan  subject  or  a structure 
and  the  horn  length  is  limited  (about  25  m)  by  cost  and  space  considerations. 

Note  also  that  a pyramidal  horn  having  a plane  floor,  walls  and  ceiling  provides 
a better  working  space  than,  for  example,  a conical  horn. 


2.2  Air  Reservoir  and  Compressor  System 

The  air  reservoir  system  for  the  Travelling-Wave  Horn  consists  essen- 
tially of  a large  cylindrical  tank  (3-4  m3)  located  inside  the  control  room 
(Fig.  11)  and  two  additional  long  cylindrical  tanks  (3-4  m3  each)  outside  the 
control  room  (Figs.  9 and  10).  The  three  tanks  are  interconnected  by  suitable 
short  pipes  of  large  diameter  to  minimize  unnecessary  pressure  differences  and 
wave  motion  between  tanks  during  the  generation  of  a simiilated  sonic  boom. 
Fiirthermore , their  combined  volume  is  sufficient  to  provide  a reasonably  constant 
reservoir  pressure  during  operation.  The  large  indoor  tank  is  connected  via  a 
short  transition  duct  directly  to  the  flap-valve  housing  (Figs.  9 and  11).  Low- 
pressvtre  air  (l  to  2 atm)  is  supplied  to  the  reservoir  by  means  of  a two-stage, 
hlgh-pressiu:e  (l  to  18  atm) , reciprocating  ccopressor  (l6  brake  horsepower, 

0.025  m3  capacity  at  85O  rpm)  and  its  small  self-contained  reservoir  (0.5  m3). 

The  high-pressure  air  in  the  con^jressor  reservoir  maintains  the  facility  reser- 
voir at  a desired,  preselected  lower  pressure  by  means  of  a sviitable  pressiire 
sensor  and  electronics,  pressure  throttle,  and  electrically  controlled, 
pneumatically  actuated  ball  valve.  Note  that  the  air  compressor  has  been  housed 
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outdoors  In  a special  acoustically  insulated  shed  (Fig.  10)  to  provide  a virtually 
noise-free  atmosphere  inside  the  control  and  test  rooms  as  well  as  the  inside  of 
the  pyramidal  horn. 

2.3  Flap  Valve  and  Control  System 

The  flap  valve  and  its  essential  requirements  for  producing  a simulated 
sonic  boom  are  illustrated  in  Fig.  15 . The  primary  function  of  the  mass-flow 
valve  (Fig.  15a)  is  to  release  in  a controlled  manner  higher  pressure  air  from 
the  reservoir  into  the  pyramidal  horn,  where  the  expanding  flow  generates  the 
simulated  sonic  boom.  To  achieve  this  simulation  each  valve  flap  must  execute 
a one-cycle,  oscillatory  rotational  motion  such  that  the  valve-throat  area  (A) 
is  approximately  a "parabolic"  function  of  time  (Fig.  15b) . For  constant 
reservoir- state  conditions  the  mass-flow  rate  (m)  of  air  through  the  valve 
throat  has  the  same  parabolic  profile  (Fig.  15b) . This  particular  mass-flow-rate 
distribution,  from  zero  to  a maximum  and  back  to  zero  again,  produces  the  desired 
N-shaped  overpressvire  signature  (/p)  of  the  simulated  sonic  boom  (Fig.  15c). 

The  mass-flow  valve,  dual  flaps  in  the  valve  housing,  three-bar  linkage, 
cam,  clutch,  flywheel,  electric  motor  and  control  system  are  shown  in  the  two 
pictures  of  Fig.  16  and  also  illustratively  in  Fig.  17.  When  the  flaps  are  in 
their  normally  closed  state  the  higher  pressure  reservoir  air  is  prevented  from 
flowing  into  the  horn.  Although  the  electric  motor  turns  the  flywheel  and  the 
nearest  part  of  the  disengaged  clutch  at  constant  speed,  the  other  part  of  the 
disengaged  clutch,  the  cam  and  valve  flaps  are  motionless.  To  actuate  the  valve 
in  order  to  generate  a simulated  sonic  boom  the  stationary  half  of  the  clutch 
is  allowed  to  move  axially  towards  the  rotating  half  and  engage  it,  causing  the 
cam  shaft  and  cam  to  suddenly  rotate  at  constant  speed.  The  rotating  cam  forces 
the  three-bar  linkage  to  synchronously  rotate  the  flaps,  first  outwards  to 
increase  the  throat  area  and  mass  flow  rate  of  air  from  the  reservoir  into  the 
horn  and  then  inwards  to  decrease  the  throat  area  and  mass  flow  rate.  After  one 
cam-shaft  revolution  the  clutch  disengages  automatically  and  also  suddenly  stops 
the  cam  and  cam-shaft  rotation.  The  flaps  have  now  returned  to  their  original 
locations  and  the  one-cycle  operation  has  ended.  The  cycle  can  be  repeated,  of 
course,  to  produce  a second  and  subsequent  simulated  sonic  booms. 

The  flap  valve  has  a centre  section  or  aerodynamic  fairing  which  has 
not  been  shown  in  Figs.  15  and  17.  A drawing  of  the  fairing  in  the  flap  valve 
is,  however,  shown  in  Fig.  18,  and  the  motion  of  the  flaps  and  various  parts  of 
the  fairing  during  the  opening  stage  of  the  valve  are  illustrated  in  Fig.  19 . 

In  Fig.  19a  the  flaps  are  shown  in  their  normally  closed  state.  The  spring- 
loaded  plates  of  the  fairing  are  pressed  together  and  a good  seal  is  maintained 
to  prevent  air  leakage  from  the  reservoir  into  the  horn.  During  the  initial 
outward  motion  of  the  flaps,  the  fairing  plates  are  forced  outwards  to  follow 
the  flaps  and  maintain  the  air  seal  (Fig.  19b).  Eventually  the  plates  are  stopped 
suddenly  by  the  fairing  front  stops,  as  the  flaps  go  beyond  the  fairing  edge 
(Fig.  19c).  At  this  time  the  transition  from  no  flow  to  flow  through  the  valve 
is  rapid,  and  it  helps  to  produce  a rapid  rise  time  on  the  front  shock  of  the 
simulated  sonic  boom.  As  the  flaps  continue  to  move  outwards,  the  first  half 
of  the  boom  is  generated  by  the  increasing  flow  of  air  into  the  horn.  The 
second  half  of  the  boom  is  generated  a,-  the  flaps  move  inwards  and  decrease  the 
flow  back  to  zero.  A rapid  transition  from  flow  to  no  flow,  when  the  flaps 
re-establish  contact  with  the  central  fairing,  helps  to  produce  a rapid  rise  time 
on  the  second  shock  of  the  boom.  Note  that  the  central  fairing  is  shaped  like  a 
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symmetrical  airfoil  (Fig.  l8)  to  minimize  flow  turbulence  and  high-frequency 
noise . 

The  two  valve  fl^s  are  each  28.0  cm  long  and  19*7  cm  wide,  and  they 
are  firmly  screwed  to  their  pivot  shafts  which  are  21.6  cm  apart  (Fig.  l8) . 

Nearly  synchronous  motion  of  the  two  flaps  is  achieved  with  a sin5)le  but  com- 
pletely adequate  three-bar  linkage  between  the  flap  pivot  shafts.  The  rotational 
motion  of  the  flaps,  which  produces  the  desired  parabolic  throat-area  variation 
and  mass-flow-rate  profile,  is  controlled  by  the  cam.  Although  the  cam  shape 
must  be  quite  precise  to  produce  the  desired  flap  motion,  the  cam  can  be  designed 
on  the  basis  of  the  predicted  shape.  When  the  flaps  are  in  the  fully  open  position 
the  valve- throat  area  can  be  as  large  as  250  cm2.  However,  the  cam  of  the  present 
system  does  not  utilize  this  whole  area,  as  the  maximum  throat  area  is  only 
156  cn^.  Note  that  the  flaps  are  quite  light  for  rapid  acceleration  response, 
and  they  have  0-ring  seals  along  their  sides  idiich  slide  over  the  valve-housing 
inner  surfaces. 

A variable- speed,  direct- current,  electric  motor  (5  brake  horsepower) 
has  the  capability  of  handling  the  average  but  not  peak  torque  requirement  of  i 

the  flap  valve.  The  reservoir  pressure  acting  on  the  flaps  during  the  valve  I 

opening  period  can  produce  a large  counter  torque.  However,  the  flywheel  ' 

(diameter  of  3I  cm,  width  of  10  cm)  adds  the  extra  capability  to  meet  the  peak 
torque  requirements.  Its  large  rotational  inertia  is  sufficient  to  give  an 
essentially  constant- speed  cam  rotation  during  the  valve  operation.  The  motor 
speed  or  the  flywheel  and  cam  speed  control  the  open-to-close  time  of  the  valve 
and  hence  the  duration  of  both  the  air  flow  through  the  valve  and  the  simulated  j 

sonic  boom.  The  flywheel  and  cam  speed  can  be  conveniently  varied  from  as  high  i 

as  600  rpm  to  100  rpm  and  even  lower  to  give  corresponding  boom  durations  as  j 

short  as  100  ms  to  500  ms  and  even  longer.  Note  that  the  reservoir  pressure  , 

dictates  the  boom  an^jlitude,  and  the  cam  shape  determines  the  boom  wave  form.  1 

An  illustration  of  the  details  of  the  fast-acting,  one-cycle,  mechanical  - 

clutch  appears  in  Fig.  20a  and  Fig.  20b.  When  the  trigger  mechanism  (Fig.  20a)  I 

is  actuated  electrically,  the  stationary  part  of  the  clutch  on  the  cam  shaft  next 
to  the  coiipling  is  released  and  forced  by  means  of  a coin5)ression  spring  along 
the  shaft  into  the  rotating  part  of  the  clutch  linked  to  the  spinning  flywheel. 

The  quick  engaging  of  the  clutch  suddenly  rotates  the  coupling  and  cam  shaft 
at  the  flywheel  speed.  As  the  cam  shaft  is  rotated  by  the  flywheel,  the  clutch 
pin  follows  the  outward  spiral  contour  of  the  clutch  housing  and  subsequently 
causes  the  clutch  to  disengage  just  prior  to  the  end  of  one  revolution  (Figs. 

20a  and  20b) . The  continuing  cam-shaft  rotation,  owing  to  rotational  inertia,  is 
then  stopped  suddenly  by  means  of  the  clutch  pin  when  it  strikes  the  rubber- 
cushioned  pin  stop.  Recoiling  rotation  of  the  cam  shaft  and  cam  is  restrained 
by  a special  spring-loaded  catch  mechanism  (Fig.  20b),  and  re-engaging  of  the 
clutch  is  prevented  by  the  trigger  mechanism  which  resets  automatically  (Fig. 

20a) . To  minimize  both  the  acceleration  and  deceleration  forces  inposed  on  the 
clutch  conponents  during  the  sudden  starting  and  stopping  phases,  the  cam,  cam 
shaft  and  associated  rotating  clutch  parts  were  designed  to  have  a small  rota- 
tional inertia. 

2.4  Reflection  Eliminator 


When  the  simulated  sonic  boom  propagates  to  the  base  of  the  horn  it 
woiild  normally  be  partially  reflected  from  the  cpen  end.  This  reflected  wave 
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would  then  disriQ)t  the  simulated  pressure  and  flow  conditions  in  the  interior 
test  section  (Fig.  9),  because  the  wave  length  of  the  simulated  sonic  boom 
can  be  as  many  as  four  times  longer  than  the  entire  horn.  In  order  to  eliminate 
or  at  least  adequately  minimize  this  imdesirable  reflection  and  its  subsequent 
echoes,  a reflection  eliminator  was  bxiilt  to  cover  the  base  of  the  horn.  The 
reflection  eliminator  is  basically  a huge  porous  piston,  as  illustrated  in 
Fig.  21.  The  porous  part  of  this  piston  is  a 2.5-cm-thick  blanket  of  microlite 
material  (12  kg/m3) . Because  the  porous  piston  is  free  to  move  on  a special 
roUer-and-track  support,  it  can  respond  to  the  incident  simulated  sonic  boom, 
being  accelerated  by  the  drag  forces  of  the  air  flowing  throiagh  the  porous 
microlite.  Significant  air  leakage  around  the  porous  piston  is  prevented  by 
enclosing  the  piston  at  its  periphery  with  a special  skirt  which  is  attached 
to  the  base  of  the  horn.  The  correct  flow  resistance  of  the  porous  microlite 
and  correct  piston  weight,  as  well  as  the  ability  of  the  piston  to  move, 
provides  a matching  of  the  impedance  of  the  duct  exit  to  that  of  the  incident 
simulated  sonic  boom,  thereby  eliminating  the  reflected  wave  (Ref.  10). 


3.  FERFORMAUCE  OF  THE  TRAVELLD^-WAVE  HORN 

The  Travelling-Wave  Horn  utilizing  the  flap  valve  has  proven  experi- 
mentally to  be  a practical  facility  for  the  simulation  of  a sonic  boom  from 
either  a SST  or  stipersonic  military  aircraft,  because  the  air  compressor  and 
reservoir  system,  flap  valve,  pyramidal  horn  and  reflection  eliminator  all 
essentially  function  correctly  as  designed.  The  simulated  sonic  boom  in  the 
horn  has  an  N- shaped  overpressure  signature  for  the  following  reasons.  Firstly, 
the  electric  motor,  flywheel,  clutch,  cam  and  three-beir  linkage  move  the  flaps 
correctly  to  give  the  desired  parabolic  throat-area  variation.  Secondly,  the 
reservoir  volume  is  siafficiently  large  to  maintain  a nearly  constant  driving 
pressiore  during  the  valve  operation  such  that  the  flow  through  the  valve  has 
the  desired  corresponding  parabolic  mass- flow-rate  variation,  which  generates 
a boom  having  the  reqiiired  N- shaped  overpressure  signature.  Thirdly,  the 
reflection  eliminator  adequately  miinimiizes  the  wave  reflection  from  the  large 
end  of  the  horn  such  that  the  N-wave  flow  and  pressxire  conditions  in  the  interior 
test  section  are  not  disripted.  Foiirthly,  jet  noise  supeiposed  on  a simulated 
sonic  boom  is  not  a problem,  because  the  large  throat  area  of  the  fl^  valve 
can  easily  pass  the  required  flow  rate  at  a low  flow  speed  and  consequently 
generates  little  jet  noise.  Note  that  more  jet  noise  than  is  desirable  is 
siperposed  on  high- amplitude  and  long-duration  booms,  which  are  only  infrequently 
needed  for  sonic-boom  tests,  but  even  this  noise  could  be  adequately  diminished 
by  utilizing  a large  cam  and  thus  a larger  valve  throat  area. 

The  N-wave  duration  requirements  are  more  than  sufficiently  satisfied 
because  the  easily  adjusted  motor  speed  can  be  set  to  give  valve  open-to-close 
times  and  thus  boom  durations  from  100  to  300  ms  and  even  longer.  Owing  to  the 
flexibility  of  setting  the  reservoir  pressure  from  as  low  as  1 atm  (absolute) 
to  as  high  as  2 atm,  the  facility  has  the  capability  of  producing  an  N-wave  having 
a lower,  equivalent  or  higher  amplitude  than  that  of  an  actual  sonic  boom.  Peak 
overpressures  as  high  as  1200,  600  and  400  N/m^  can  be  achieved  easily  for 
respective  N-wave  dtirations  of  100,  200  and  300  ms.  The  rise  times  of  the  front 
and  rear  shocks  of  the  simulated  sonic  boom  are  not  constant  but  depend  very 
weakly  on  N-wave  anplitude  and  more  strongly  on  duration.  For  exanple,  a typical 
rise  time  of  the  front  shock,  defined  as  1.25  times  the  time  for  the  overpressure 
to  rise  from  10  to  90^  of  its  peak  value,  is  about  6,  8 and  10  ms  for  corresponding 


7 


i 

I 

i 

[ 

i 


ruations  of  100,  200  and  3OO  ms.  The  rise  time  of  the  rear  shock  is  always 
somewhat  longer  (lO  to  2%). 


The  Travelling-Wave  Horn  has  beer,  designed  for  continuous  operation 
such  that  a sufficiently  rapid  sequence  of  simulated  sonic  booms  can  b e produced 
to  facilitate  structural  response,  fatigue  and  damage  studies.  The  number  of 
booms  per  unit  time  depends  on  both  the  mass  of  air  discharged  from  the  reservoir 
per  boom  and  the  capacity  of  the  con^jressor  to  refill  the  reservoir.  For  the 
case  of  a long-duration  N-wave  of  300  ms  having  a modest  peak  oveipressure  of 
100  n/u^,  the  generation  rate  has  been  found  to  be  about  three  booms  per  minute. 
For  a shorter  duration  boom  of  200  ms  having  the  same  aii5)litude,  the  facility 
can  easily  produce  six  booms  per  minute.  If  the  duration  is  still  shorter  at 
100  ms  the  generation  rate  increases  to  about  fifteen  booms  per  minute.  When 
the  N-wave  aii5)litude  is  doubled  from  100  to  200  N/m^,  then  the  generation  rate 
is  reduced  by  a factor  of  two. 


As  already  mentioned,  the  valve  flaps  move  correctly  and  give  a para- 
bolic throat-area  variation  with  time.  Although  the  area  variation  cannot  be 
measured  directly  diiring  the  generation  of  a simulated  sonic  boom  to  substantiate 
the  preceding  statement,  it  can  be  inferred  directly  from  a measurement  of  the 
angular-displacement  history  of  a flap.  Note  that  the  throat  area  is  directly 
proportional  to  the  angular  displacement  of  f flap,  provided  that  the  rotation  is 
sufficiently  small.  Measured  angular-displacement  histories  of  the  top  flap, 
which  are  virtually  the  same  as  those  for  the  bottom  flap,  are  shown  in  Fig.  22. 
Results  are  given  not  only  for  three  different  N-wave  durations  of  100,  200  and 
300  ms,  but  also  for  three  different  reservoir  overpressures  of  6.9,  13.8  and 
20.7  kN/m2.  These  measured  profiles  have  been  reproduced  in  Fig.  23  where  they 
can  be  better  ccjnpared  to  a parabolic  profile.  From  these  and  other  angtilar- 
displacement  histories  it  can  be  concluded  that  the  cam  shape  and  associated 
flap  rotation  produce  the  desired  parabolic  throat-area  variation. 

The  capability  of  the  Travelling-Wave  Horn  for  generating  both 
different  aii5)litude  and  duration  simulated  sonic  booms  having  a good  N- shaped 
signatiire  is  aptly  demonstrated  by  the  measiu*ed  overpressure  signatures  shown 
in  Figs.  24,  25  and  26.  These  N- shaped  signattires  are  qxute  representative  of 
those  of  actual  sonic  booms. 


The  measiu’ed  signatures  also  show  that  low-amplitude  and  short-duration 
booms  are  virtually  free  of  jet  noise,  and  significant  jet  noise  only  occurs 
siperposed  on  infrequently  used  booms  having  a high  amplitude  and  long  duration. 
It  is  now  worthwhile  to  conpare  the  virtually  noise-free  signatures  of  Figs.  24, 
25  and  26  for  the  case  of  the  flap  valve  to  the  noisy  signatures  of  Figs.  7 and 
8 for  the  other  case  of  a plug  valve  and  pyramidal  horn  with  and  without  a jet- 
noise  absorber.  From  this  conparison  it  is  qvdte  obvious  that  the  facility 
modification  from  a plug  valve  to  a much  larger  flap  valve  was  more  effective  in 
reducing  the  jet  noise  than  the  addition  of  a jet-noise  absorber  to  the  original 
facility  utilising  a plug  valve. 

Some  additional  high-frequency  noise  appears  ahead  of  the  front  shock 
of  short- duration  simulated  sonic  booms,  as  can  be  seen  in  Fig.  24.  Such 
undesirable  precursor  noise  is  due  to  mechanical  banging  of  mainly  clutch  parts 
during  the  fast  clutch  engaging  process.  This  noise  is  of  course  more  severe 
for  higher  flywheel  and  cam-shaft  speeds  required  for  short- duration  N-waves. 
Precursor  noise  is  not  much  of  a problem  because  it  is  virtually  absent  from 
overpressure  signat\ires  having  a dviration  greater  than  about  100  ms  (see  Figs. 

25  and  26). 


I 


1 


8 


The  low-frequency  perturbations  in  the  basic  N- shaped  overpressure 
signatures  shown  in  Figs.  24,  25  and  26  are  not  due  to  a defective  reflection 
eliminator,  but  result  instead  from  enclosing  the  large  end  of  the  horn  and 
reflection  eliminator  in  the  test  room  (Fig.  9) • The  simulated  sordc  boom,  on 
passing  out  of  the  horn  through  the  porous  piston  into  the  test  room,  cause 
low-frequency  wave  motion  in  the  test  room  (Helmholtz  resonator  effect) . Such 
waves  siibsequently  enter  the  horn  through  the  porous  piston  and  slightly  disrupt 
the  flow  and  pressure  conditions  in  the  interior  test  section,  thereby  causing 
undesirable  perturbations  in  the  N-shaped  signature.  To  minimize  the  enclosing 
effect  of  the  test  room,  two  large  test-room  doors  directly  behind  the  reflection 
eliminator  are  opened  diiring  normal  facility  operation. 

Before  discussing  the  rise  times  of  the  front  and  rear  shocks  of  the 
simulated  sonic  boom  it  is  in^ortant  to  first  define  the  rist  time.  The  rise 
time  is  defined  herein  to  be  1.25  times  the  time  for  the  overpressure  to  rise 
from  10  to  90^  of  its  peak  value.  The  peak  value  should  be  that  of  the  basic 
N-shaped  signature  and  thus  exclude  any  spurious  peak  caused  by  the  jet  noise 
or  perturbation  caused  by  the  test-room  enclosing  effect.  This  rise- time 
definition  gave  consistent  and  meaningful  results.  The  measured  rise  times  for 
repeated  trials  were  very  similar  because  the  10  and  90^  points  of  the  shock 
profile  were  well  defined  in  each  case  and  also  consistent,  as  was  the  rapid 
overpressure  rise  between  these  two  points.  Note  that  the  0 and  100^  points 
of  the  shock  profile  were  not  well  defined  and  varied  substantially  from  one  test 
to  the  next.  The  rise- time  definition  gave  meaningful  results  not  only  because 
the  measured  rise  times  were  cu;isistent,  but  also  because  the  subjective  loudness 
of  a shock  is  mainly  a function  of  the  rapid  rate  of  change  of  overpressure  and 
only  weakly  dependent  on  the  total  overpressure  change.  Furthermore,  in  the  case 
of  a shock  profile  which  rises  asynptotically  to  its  peak  value  (rear  shock  of  an 
N-wave) , a rise-time  definition  based  on  the  overpressure  change  from  a minimum 
(0^  point)  to  a maximum  (ICX)^  point)  would  yield  an  infinite  rise  time,  whereas 
the  definition  adopted  in  this  work  would  give  a realistic  finite  value. 

Rise  times  of  the  front  and  rear  shocks  of  simulated  sonic  booms 
shown  in  Figs.  24,  25  and  26  are  not  all  equal  or  constant,  but  depend  mainly 
on  the  valve  open-to- close  time  or  N-wave  duration  and  to  some  extent  on 
reservoir  pressure  or  N-wave  arplitude.  Neglecting  the  small  effect  of  reser- 
voir pressure,  the  rise  time  of  the  front  shock  averages  about  6,  8 and  10  ms 
for  a duration  of  100,  200  and  300  ms  respectively.  The  average  rise  time  of  the 
rear  shock  is  generally  10  to  25^  longer.  These  rise  times  of  simulated  sonic 
booms  are  certainly  longer  than  most  of  those  taken  from  measurements  of  actual 
sonic  booms.  For  example,  the  rise  time  of  the  front  shock  of  a sonic  boom  from 
the  Anglo-French  Concorde  SST  ranges  from  as  low  as  0.1  ms  to  5 &nd  sometimes 
even  longer  (Ref.  23),  and  probably  has  an  average  of  about  1 ms.  On  the  other 
hand,  the  average  rise  time  of  638  sonic  booms  from  militaxy  aircraft  was  found 
to  be  4 to  5 ms  (Ref.  24). 

The  rise  time  of  the  simulated  sonic  boom  from  the  Travelling-Wave 
Horn  is  certainly  short  enough  at  6 to  10  ms  to  facilitate  most  room-resonance 
studies  and  structural-panel  response  and  damage  investigations.  However,  most 
human-response  studies  (startle,  hearing,  heart  rate,  sleep)  and  also  animal- 
response  tests  (cochlea  damage)  require  a sonic  boom  having  a much  shorter  rise 
time.  Additionally,  it  is  inportant  to  control  the  rise  time  and  be  able  to 
vary  it  from  1 ms  or  even  shorter  if  possible  to  10  ms  and  longer.  Although 
the  ccoplementatry  ITEIAS  Loudspeaker-Driven  Booth  can  generate  a simulated  sonic 
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boom  having  a rise  time  which  can  be  controlled  from  1 to  10  ms  and  even  longer, 
in  order  to  facilitate  htonan  and  animal -response  studies,  it  was  hc^ed  that  the 
UTIAS  Travelling- Wave  Horn  wovild  have  a similar  capability. 

Much  thou^t  was  originally  given  to  the  design  of  the  flap  valve  to 
achieve  a short  rise  time  on  both  shocks  of  the  simulated  sonic  boom,  hopefully 
as  short  as  1 ms.  Furthermore,  a sin5)le  scheme  was  contrived  to  control  the 
rise  time,  such  that  it  could  be  varied  from  1 to  10  ms  or  longer  as  desired. 

The  first  central-fairing  design,  shown  schematically  in  Fig.  27,  feattired  two 
very  thin  steel  plates.  As  the  valve  flaps  initially  move  outwards  (Fig.  27a) 
and  break  free  of  the  central  fairing,  the  steel  plates  are  simultaneously 
translated  and  rotated  otitwards  from  their  normally  retracted  position  into 
their  outermost  location  (Fig.  27b).  Later,  when  each  flap  retvims,  each 
associated  steel  plate  is  quickly  retracted  to  clear  the  space  for  the  flap  and 
thus  avoid  severe  flap  iir^jact  with  the  plate.  In  experiments  with  this  first 
central  fairing,  short  rise  times  of  about  1 ms  on  the  front  shock  and  2 to  3 
ms  on  the  rear  shock  were  achieved,  provided  that  the  steel  plates  moved  very 
quickly  into  their  outermost  location  and  later  back  into  their  retracted  position. 
On  the  other  hand,  when  the  plates  were  not  allowed  to  move  and  fixed  in  their 
retracted  state,  both  rise  times  were  very  long  at  about  15  ms.  By  siaply 
controlling  how  qiiickly  the  steel  plates  could  be  extended  or  retracted,  it  was 
fovind  that  the  rise  time  of  the  front  shock  could  be  correspondingly  controlled 
and  varied  from  1 to  15  ms.  Also,  the  rise  time  was  not  strongly  affected  by  the 
N-wave  duration  or  anplitude.  These  experiments  were  very  inportant  in  that  they 
confirmed  the  two-fold  concept  of  generating  a simulated  sonic  boom  having  a 
short  rise  time  of  about  1 ms  and  controlling  the  rise  time  with  a special  control- 
fairing  mechanism. 

Difficulties  were,  however,  encountered  with  the  first  control- fed  ring 
design  (Fig.  27) , which  unfortunately  prohibited  its  use  for  normal  facility 
operation.  On  the  return  motion  of  the  flaps  to  the  central  fairing,  the  flaps 
would  occasionally  strike  the  steel  plates  with  sufficient  inpact  to  cause 
permanent  plate  bending.  This  bending  would  often  cause  other  valve  parts  to 
bind  and  sometimes  fail.  Such  valve  failure  occurred  as  frequently  as  once  for 
every  20  to  50  simulated  sonic  booms,  and  valve  repair  generally  required  a few 
man  days  of  effort.  To  overcome  this  problem  a second  central  fairing  (Figs. 
l8  and  19),  which  was  used  for  all  work  given  in  this  report,  was  designed  and 
tested.  This  second  central  fairing  has  never  failed  because  each  of  the  two 
plates  is  firstly  hinged  at  the  centre  to  reduce  inpact  forces  and  secondly 
thicker  to  avoid  bending  failure.  However,  the  penalty  of  using  this  second 
central  fairing  has  been  that  the  rise  time  of  the  front  and  rear  shocks  of  the 
simulated  sonic  boom  are  xaidesirably  long  at  6 to  10  ms  and  longer.  The  thick, 
hinged  plates  of  the  central  fairing  vinfortunately  do  not  move  sufficiently  fast 
or  in  just  the  correct  manner  to  produce  the  rapid  starting  and  stepping  flow 
processes  needed  to  produce  short  rise  times  on  the  front  and  rear  shocks 
respectively.  If  simulated  sonic  booms  having  a short  rise  time  of  less  than 
6 ms  are  required  for  normal  facility  operation,  then  modification  or  possibly 
new  design  of  the  central  fairing  and  valve  flaps  is  necessary. 

AJIALYSIS  FOR  THE  TRAVELLIHG-WAVE  HORN 


A combination  of  gasdynamlc  and  acoustic  theory  is  used  to  successfully  ' 

predict  the  main  featvires  of  the  simulated  sonic  boom  (wave  form,  amplitude, 


duration,  rise  time).  Gasdynamic  theory  is  used  first  to  describe  the  time- 
varying  conditions  of  the  air  in  the  reservoir  (e.g.,  pressure  and  ten5)erature) 
and  the  mass  flow  rate  of  air  out  of  the  reservoir  through  the  opening  and 
closing  flap  valve  into  the  pyramidal  horn.  Acoustic  theory  is  then  en^ployed 
to  describe  the  resulting  wave  motion  in  the  horn.  Because  of  certain  assuii5)- 
tions  made  in  the  gasdynamic  analysis,  the  mass-flow  rate  of  air  out  of  the 
reservoir  can  be  determined  without  solving  for  the  wave  motion  ir.  the  horn, 
and  it  thus  serves  as  a convenient  boundary  condition  for  the  st.parate  acoustic 
analysis.  Although  this  type  of  analysis  for  a travelling-wave  horn  is  based 
on  work  given  originally  in  Refs.  6 and  7 and  much  more  con5)letely  in  Ref.  10, 
the  present  analysis  features  not  only  a new  interesting  method  of  solution  but 
also  new  results. 

^.1  Gasdynamic  Analysis  of  the  Reservoir 

The  adr  reservoir  of  the  Travelling-Wave  Horn  has  a finite  volxane 
(10.2  m3).  Consequently  the  properties  of  the  reservoir  air  are  not  fixed  with 
time  but  vary  continuously  as  the  air  discharges  out  through  the  flap  valve  and 
generates  the  simulated  sonic  boom  in  the  pyramidal  horn.  In  order  to  predict 
the  reservoir  conditions,  assume  that  the  state  properties  axe  spatially  uniform 
and  all  changes  with  time  are  adiabatic  and  reversible.  Then  the  following 
isentropic  expressions  can  be  used  to  relate  the  time- dependent  reservoir  tenper- 
at\ire  To(t),  sound  speed  ao(t) , pressure  po(t)  and  density  po(t)  to  their  initial 
fixed  values  which  are  denoted  by  an  overhead  bar. 
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The  symbol  y denotes  the  ratio  of  the  specific  heat  at  constant  pressure  (Cp) 
to  that  at  constant  volume  (CV) • Now  assume  that  the  air  flow  from  the  reservoir 
through  the  flap  valve,  which  is  essentially  a convergent-divergent  duct  flow 
with  a slowly  varying  throat  area,  is  not  only  adiabatic  and  reversible  but 
also  one-dimensional  and  quasisteady.  Isentropic  expressions  given  below  then 
relate  the  time-dependent  properties  of  the  reservoir  air  to  those  at  the  valve 
throat  (denoted  by  the  subsciprt  *)  and  also  to  the  flow  Mach  number  M^(t). 
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The  following  equation  for  the  mass-flow  rate  of  air  per  unit  area  at  the  valve 
throat  m».(t)  can  be  expressed  as  a function  of  the  Mach  nunber  M^(t) , which  is 
sirply  the  ratio  of  the  flow  speed  u^(t)  and  sound  speed  a*(t) . 


m^(t)  - p*(t)  a^(t)  M^(t) 


(4) 


If  the_flow-area  variation  A^(t)  at  the  valve  throat  is  expressed  as  A^N^(t), 
where  A^  is  the  maximum  throat  area  and  N^(t)  describes  the  area  variation, 
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then  toe  mass  flow  rate  is  m^(t)  A^N^(t)  or  p^(t)  a„(t)  JL(t)  A^N^(t) . For  this 
mass  flow  rate  of  air  out  of  the  reservoir  having  a volume  V,  the  resulting  rate 
of  decrease  of  reservoir-air  mass  is  sin^ply  Vp'(t),  where  the  prime  denotes 
‘^^^®^®iitiation  with  respect  to  time.  The  reservoir  density  is  thus  governed 
by  the  following  differential  equation. 


p;(t) 


P*(t)  a^(t)  I^(t)  Ag^N^(t)/V  = 0 


(5) 


In  order  to  solve  this  differential  equation  in  conjunction  with  previous 
expressions,  toe  equation  of  state  (p  = pET)  and  sound-speed  expression 
(a  = 7RT),  a boundary  condition  for  toe  flow  at  the  valve  throat  is  needed. 

first  the  case  of  choked  flow  at  the  valve  throat  for  which 
I^(t)  does  not  change  with  time  but  equals  unity.  This  is  toe  required  boundary 
condition  if  the  flow  is  choked.  Choked  flow  can  be  assumed  to  persist  with 
time  as  long  as  the  throat  pressure  P^(t)  exceeds  the  ambient  pressure  in  the 
pyramidal  horn,  toich  can  be  taken  to  be  atmospheric  pressure  (p)  . Or,  choked 
flow  occurs  when  toe  reservoir  to  atmospheric  pressure  ratio  Po(t)/p  exceeds 
the  reservoir  to  throat  pressure  ratio  Po(t)/p„(t).  From  Eq.  3 with  M^(t) 
equal  to  unity,  this  critical  pressure  ratio  Po(t)/p„(t)  is  constant  and 

equal  to  (7+l)/2  raised  to  the  power  of  For  the  case  of  air 

with  7 equal  to  1.4  the  critical  pressiire  ratio  is  thus  I.892  or  toe  reciprocal 
01  U • 

By  Using  Eqs.  2 and  3 letting  I^(t)  equal  ixnity,  the  differential 
equation  (Eq.  5)  can  be  expressed  in  terms  of  only  one  unknown  TgCt),  ao(t), 
Po(t)  or  po(t)  and  integrated  e2q)licitly.  The  sirplest  procedure  is  to  use 
the  unknown  ao(t)  to  obtain  the  following  result. 
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Other  results  for  To(t),  Po(t)  and  po(t)  follow  from  Eq.  2.  The  mass  flow  rate 
per  mit  area  ii^(t),  normalized  by  the  product  of  atmospheric  density  p and  sound 
speed  a,  can  be  derived  easily  from  the  preceding  resialts,  and  it  is  given  below. 
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Note  that  the  nondimensional  parameter  ^^to/V  and  the  integral  of  N (t)  deter- 
mine how  fast  the  reservoir  conditions  and  mass  flow  rate  per  unit  area  change 
with  time.  The  changes  are  rapid  and  large  if  N„(t)  is  unity  and  the  parameter 
is  large,  corresponding  to  a small  air  reservoir  (V),  large  valve  (A„)  and  a 
long  valve  open-to-close  time  (to) . * 


Now  consider  the  other  case  of  unchoked  or  subsonic  flow  at  the  valve 
throat  for  which  M^( t)  is  less  than  unity  and  not  constant.  For  this  case  a 
reasonable  boundary  condition  for  the  flow  at  the  valve  throat  is  to  take  the 
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static  pressvtre  p^(t)  as  constant  and  eqvial  to  atmospheric  pressure  p.  An 
Immediate  consequence  which  follows  from  Eqs.  2 and  3 is  that  all  other  throat 
conditions  such  as  T^(t),  a*(t)  and  p*(t)  are  also  constant.  By  xising  this 
assim?)tion  and  Eqs.  2 and  3>  the  previous  differential  equation  for  the  density 
(Eq.  5)  can  now  be  ejcpressed  in  tenns  of  only  one  unknown  M^(t)  as  shown  below. 
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Once  l^(t)  is  obtained  from  a solution  to  this  differential  equation  the 
reservoir  time-dependent  conditions  and  mass  flow  rate  follow  from  Eqs.  2,  3 
and  4. 


Eqviation  8 can  be  solved  exactly  for  only  certain  values  of  y (e.g., 
5/3 j 3/2  and  7/5).  For  the  relevant  case  of  air  (7  = 7/5) » the  exact  solution 
can  be  given  as  follows, 
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where  x equals  ^L(t)/•^  and  the  initial  value  x is  The  initial  flow 

Mach  number  at  time  t equal  to  zero  follows  from  Eq.  3 with  p^(t)  replaced 
by  p as  shown  below. 


Note  that  a^  in  Eq.  9 ts  the  sound  speed  of  the  air  at  the  valve  throat  and  not 
in  the  reservoir.  However,  it  is  related  to  the  initial  values  of  reservoir 
soimd  speed  and  pressure  by  Eqs.  2 and  3,  as  shown  below. 


(11) 


Now,  for  a given  initial  reservoir  sound  speed  a-Q,  pressure  ratio  Po/p  and 
volume  V,  valve  open-to-close  time  or  flow  d\iration  t©,  and  throat-area 
distribution  A^N^(t),  the  flow  Mach  nuntoer  M^(t)  can  be  obtedned  from  Eq.  9 
by  an  iterative  procedure. 


The  noimalized  ftmction  N„(t)  for  the  throat-area  variation  has  thus 
far  been  left  in  a general  form  suc5  that  any  actual  or  idealized  form  of  N^(t) 
can  be  used  in  the  analysis.  At  this  stage,  however,  it  is  worthwhile  to 
introduce  an  ideal  "parabolic"  form  for  N^(t),  which  is  given  below. 
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N„(t)  = 


Hence,  the  integral  in  Eq.  9 is  given  as  follows, 


t < 0 


0 < t < t. 


t < t 

O 


/ W-|[3(l;)  -<f) 


and  it  has  a maxlimm  value  of  2/3.  It  should  be  noted  that  in  practice  the 
cam  shape  and  resulting  valve-flap  motion  are  designed  to  produce  this  parabolic 
throat-area  variation  (Eq.  12),  in  order  to  generate  a simulated  sonic  boom 
having  a duration  to* 

Once  M,(t)  has  been  determined  by  using  Eqs.  9,  10,  11  and  12  or  13  for 
a fixed  initial  reservoir  pressxire  Pq  and  sound  speed  ao>  the  reservoir  conditions 
with  time  follow  from  Eqs.  2 and  3*  For  example,  the  expression  for  the  reservoir 
pressure  is  given  below. 


Po(‘) 


l^(t) 


Furthermore,  by  using  Eq.  3 and  previous  results,  the  mass  flow  rate  per  unit 
area  m^(t)  can  be  expressed  in  the  following  nondimensi onal  form. 


a r Fq  "^7 

f i^(t) 

pa  a P 


Although  it  may  not  be  readily  apparent  from  Eqs.  9>  and  15,  the  nondimensional 
parameter  a^A^to/V  in  the  solution  for  J^(t),  Po(t)  and  n^(t)  determines  how  fast 
the  reservoir  conditions  and  mass  flow  rate  change  with  time,  like  the  previous 
case  for  choked  flow  at  the  valve  throat. 

In  past  work  with  the  Travelling-Wave  Horn  utilizing  a plug  valve 
(Ref.  10)  the_reservoir  pressure  was  most  often  sufficiently  high  to  produce 
choked  flow  (po  > 1.89p).  In  some  cases,  however,  the  reservoir  pressure  would 
decrease  sufficiently  such  that  the  flow  would  become  unchoked  at  later  times. 
When_the  reservoir  pressure  falls  to  1.89p,  the  new  reservoir  conditions  (fo> 
ao>  Po  and  po)  at  this  time  can  be  determined  from  the  choked- flow  solution  and 
then  used  as  input  for  the  unchoked- flow  solution.  This  procedvure  may  be  tedious 
but  it  is  not  laadvily  difficult.  In  the  present  work  for  a large  flap_valve, 
however,  the  reservoir  pressure  is  always  sufficiently  low  (p^  < 1.89p)  such 
that  the  two  cases  of  choked  flow  and  a combination  of  choked  and  unchoked  flow 
are  of  no  concern  to  this  work.  Henceforth  the  work  will  pertain  solely  to  the 
case  of  unchoked  flow  at  the  valve  throat. 

For  the  case  of  unchoked  flow  it  is  inconvenient  in  the  analysis  to 
obtain  ^^(t)  iteratively  from  the  exact  solution  (Eq.  9)  For  7 equal  to  1.4,  or 


r..i 


numerically  from  the  original  differential  equation  (Eq.  8)  for  other  values  of  7. 
Fortunately,  this  inconvenience  can  be  circumvented  by  using  an  approximate  but 
very  accurate  ea^licit  expression  for  \(t).  This  approximate  expression  inposes 
no  restriction  on  the  value  of  7 for  normal  values  between  1 and  5/3.  To  obtain 
this  expression  the  variable  coefficient  in  Eq.  8,  namely  1 + /2  raised 

to  the  power  of  (2-7)/(7-l),  is  first  expanded  in  a binomial  series.  Such  an 
expansion  is  valid  because  (7-l)^^^(t)/2  is  a small  parameter  less  than  unity 
which  never  exceeds  (7-l)/2  (e.g.,  0.2  for  air).  Then  the  differential  equation 
can  be  integrated  easily  to  give  the  following  result. 


M^(t)  - \ [^^(t)  - 5^]  + ~ [^^(t) 


+ • • • + 


u p u 

-ir-^  J = 0 

o 

Now  let  an  explicit  solution  for  M^(t)  take  the  following  form, 

« 00 

^(t)  = % ^ a.aV(t), 


(16) 


(17) 


i=0 


a = 

I(t)  =1  J N*(y)dy, 


vdiich  is  reasonable  provided  that  ai(t)  is  a small  parameter  having  an  absolute 
value  less  than  unity.  Note  that  in  the  present  work  I(t)  never  exceeds  unity 
and  a typically  varies  from  0 to  0.2. 


When  Eq.  17  for  >L(t)  is  substituted  into  Eq.  16  and  all  terms  of  third 
ordei  and  lower  are  retained,  the  following  approximate  results  for  M„(t)  can  be 
obtained.  ^ 


M»(t)  _ _ . _ 

= 1 + a^ai(t)  + a^o:  I (t)  + a2a-^I-^(t) 


(18) 


tA, 


a 


= i(t)  = I J ° %(y)dy 

o 

^ -4 

^ 3%[2  + (2  - 7)wg] 
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a -16(2  - y) 

2 9\[2  + (2  - 

64(2  - y)[2  - 5(2  - y)^] 

2 ■ 8B^[2  + (2  - y)i^]9 

This  ejqjlicit  exjjression  for  M^(t)  can,  of  course,  be  used  in  Eqs.  l4  and  15  to 
evaluate  the  reservoir  pressure  Po(t)  and  mass  flow  rate  per  \anit  area  m^(t). 

Results  of  the  unchoked- flow  analysis  are  illustrated  graphically  in 
Figs.  28  to  33*  In  each  case  the  throat-area  variat,'’  n with  time,  N^(t),  was 
taken  to  be  parabolic  (Eqs.  12  and  13).  Consider  fii„t  the  results  for  the  flow 
Mach  nianber  Mx.(t)  given  by  Eq.  iB  and  shown  in  Fig.  28.  Three _plots  for  three 
different  but  typical  veilues  of  the  nondimensional  grouping  a A„to/V  are  shown 
in  this  figure.  In  each  plot  four  different  variations  of  M^tt)  appear,  corre- 
spond!^ to  four  different  initial  values  of  the  reservoir  to  atmospheric  pressure 
ratio  Po/p  (or  equivalently  the  Mach  number  M,^) . From  these  results  it  can  be 
seen  that  J^(t)  changes  from  its  initial  value  more  markedly  when  Po/p  (or 
is  lower  and  also  when  the  value  of  a^A^tQ/v  is  larger  (large  throat  area, 
long  valve  open-to-close  time,  small  reservoir  volume).  It  should  be  noted  here 
that,  for  normal  operation  of  the  Travelling- Wave  Horn  to  generate  a simulated 
sonic  boom  having  an  anjjlitude  from  20  to  200  N/m^  and  a duration  from  100  to  300 
ms,  the  change  in  ^^(t)  is  small  (less  than  20^)  because  the  reservoir  was  inten- 
tionally designed  to  have  a large  volume. 

It  is  worth  mentioning  that  the  flow  Mach  number  J^(t)  can  diminish  to 
zero  if  the  valve  open-to-close  time  is  sufficiently  long.  One  such  result  is 
shown  for  interest  in  Fig.  28c  (pq/p  = 1.007)  . Once  the  Mach  nuniber  diminishes 
to  zero  the  flow  from  the  reservoir  ceases,  even  though  the  valve  may  still  be 
partly  open  and  Eq.  l8  predicts  a negative  Mach  number  or  reverse  flow.  It 
shovild  be  noted  that  there  is  no  mechanism  in  the  present  quasisteady  gasdynamic 
analysis  to  predict  a reverse  flow  or  flow  oscillations  like  that  of  a second- 
order  differential  equation  describing  a Helmholtz  resonator  (spring,  mass  and 
dar5)er  system)  . 

The  mass  flow  rate  per  unit  area  at  the  valve  throat,  ii^(t),  changes 
with  time  like  the  Mach  number  l^(t)  shown  in  Fig.  28,  because  it  is  directly 
proportional  to  M)(.(t)  (see  Eq.  15).  Consequently,  these  results  are  not  shown 
graphically.  Of  more  interest,  however,  is  the  variation  of  the  mass  flow  rate 
™*(t)A^N^(t) . Although  N^(t)  is  a parabolic  function  of  time  (Eq.  12),  the  mass 
flow  rate  will  deviate  from  this  ideal  profile  because  Ji^(t)  is  not  constant  but 
varies  like  l^(t) . The  degree  of  deviation  for  different  values  of  Pq/p  and 

to/y  can  readily  be  seen  in  Fig.  29,  where  the  nondimensional  mass  flow  rate 
™,(.(t)/pa  is  plotted  versus  time  t/tg  and  compared  directly  to  the  corresponding 
parabolic  profile  given  by  m^N^(t)/^  [El  is_the  initial  value  of  mj(^(t)].  For  a 
small  value  of  po/p  and  a large  value  of  a^A^to/V,  it  can  be  seen  that  the  mass 
flow  rate  deviates  more  markedly  from  the  ideal  parbolic  profile. 

The  variation  of  the  nondimensional  pressure  Po(t)/p,  density  Po(t  )/p 
and  ter5)erabure  To(t)/T  with  nondimensional  time  t/to,Tihich  are  each  given  by  a 
combination  of  Eqs.  2,  l4  and  l8,  are  shown_in_Figs._39,  31  and  32  respectively. 
The  marked  effects  of  large  values  of  both  Po/p  and  a^A^tQ/v  on  making  the 
variation  large  can  readily  be  seen.  However,  by  noting  that  the  ordinate  scale 
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has  been  greatly  expanded,  it  can  also  be  seen  that  the  total  relative  change 
in  each  Po(t)/p,  Po(t)/p  and  T(t)/T  profile  is  normally  small,  even  for  the 
largest  values  of  Pq/p  and  a^A^to/V. 

The  results  for  the  reservoir  pressure,  density  and  ten^jerature  shown 
in  Figs.  30 > 31  and  32  do  not  bring  out  certain  information  pertinent  to  the 
selection  of  the  reservoir  volume  for  a sonic-boom  simulator.  Since  the  driving 
force  for  the  air  flow  from  the  reservoir  into  the  horn  is  not  directly  due  to 
the  absol\ite  reservoir  pressvire  Po(t)  or  pressure  ratio  Po(t)/p  but  rather  to 
the  pressure  difference  PoCb)  - p,  it  is  iinportant  that  this  pressure  difference 
remains  essentially  constant  during  the  generation  of  a simulated  sonic  boom. 
Results  for  the  nondimensional  pressure  difference  (po(t)  - p) /(Po  **  P)  versus 
time  t/to  are  shown  in  Fig.  33*  They  exhibit  the  same  temporal  variation  and 
dependence  on  a^A^to/v  as  the  previous  resiilts  for  Po(t)/p  (Fig.  30)  but  an 
opposite  dependence  on  Po/P-  The  change  in  (poCt)  - p)/(Po  - P)  with  time  is 
larger  for  smaller  values  of  Pq/p*  Consequently,  the  reservoir  volume  for  a 
facility  must  be  sufficiently  large  or  a^A^to/V  sufficiently  small  such  that 
(Po(t)  - p)/(Po  - p)  varies  insignificantly  for  the  lowest  value  of  the  reser- 
voir pressure  Po  which  one  expects  to  use. 

4 .2  Acoxistic  Analysis  of  the  Simulated  Sonic  Boom 

The  sound  or  simulated  sonic  boom  in  the  pyramidal  horn  is  generated 
by  the  influx  of  air  at  the  horn  apex.  The  well-known  spherical  wave  equation 
which  is  given  below  is  enployed  to  describe  the  wave  motion  in  the  horn  (Ref.  25) 

(19) 

at  dr 

The  respective  symbols  (p  and  r denote  the  velocity  potential  and  the  radial 
distance  measured  from  the  horn  apex.  Becaxise  only  an  outward  moving  wave  is 
generated  by  the  air  entering  the  horn  and  the  reflection  of  this  wave  at  the 
large  end  is  neglected  or  considered  eliminated  by  the  reflection  eliminator, 
the  general  solution  of  the  wave  equation  can  be  expressed  in  the  following 
simple  form. 


4,  = f(T)/r 


(20) 


T = t.-  (r  - r^)/a 


The  symbol  Tq  denotes  the  radial  distance  (I.58  m)  at  which  the  flap  valve  is 
attached  to  the  pyramidal  horn,  and  t is  a retarded  time  having  an  initial  value 
of  zero  at  the  wave  front. 

The  overpressure  ^ and  particle  velocity  ^ of  the  wave  in  the  horn 
are  related  to  the  velocity  potential  4,  and  also  to  the  general  function  f(T) 
which  describes  the  wave  signature,  as  illustrated  below. 

^ =-  = (21) 
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Au  = 


(22) 


- f'(T)  - fir) 

a r r 


I 
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The  prime  denotes  differentiation  of  the  variable  with  respect  to  the  argument 
given  in  the  following  brackets.  It  is  worth  noting  that  the  overpressure  is 
directly  proportional  to  the  derivative  of  f(T)  and  decays  with  distance  like 
l/r.  The  particle  velocity,  however,  is  essentially  directly  proportional  to 
f(T)  at  small  radii  or  in  the  so-called  near- field  and  diminishes  like  l/r^, 
whereas  at  Istrge  radii  or  in  the  far-field  it  is  directly  proportional  to 
f'(T)  and  decays  like  l/r.  Hence,  the  overpressvire  and  particle- velocity 
signatiires  have  the  same  shape  only  in  the  far-field. 

In  order  to  determine  the  function  f( t)  and  thus  the  overpressure 
and  particle  velocity  of  the  wave  in  the  horn,  the  mass  flow  rate  at  the  valve 
throat  from  the  gasdynamic  analysis,  W^(t)A  N^(t) , is  used  as  the  boundary 
condition  for  the  acoustic  analysis.  Firstly,  it  is  assumed  that  the  flow  is 
quasisteady  from  the  valve  throat  to  the  downstream  location  where  the  flap 
valve  is  joined  to  the  horn.  This  location  is  at  a radial  distance  ro,  and 
the  corresponding  area  is  denoted  by  Aq.  The  mass  flow  rate  at  this  location 
is  or  pAqAu  to  first  order.  Then,  this  mass  flow  rate  ceui  be  equated 

to  nL(t)A^N^(T)  from  the  gasdynamic  analysis,  and  An  can  be  eliminated  from  the 
resulting  expression  by  using  Eq.  22,  to  give  the  following  first-order,  linear 
differential  eqviation  for  f(T). 

- ar  A 

f*(T)  f(T)  = - — ^ (23) 

O P 

When  the  mass  flow  rate  ni^(T)A^N^(  t)  is  zero  before  time  r equal  to  zero,  the 
fiaiction  f(T)  and  its  derivatives  are  all  equal  to  zero.  For  a mass  flow  rate 
starting  at  t equal  to  zero  and  ending  at  t equal  to  to>  the  solutions  for  f(T) 
and  f'(T)  during  this  time  interval  can  be  obtained  easily  from  Eq.  23,  and  they 
are  given  below. 


ar  A / “ \ / - \ 

f(T-)  ^ — exp  ^ jJ  n^(y)N^(y)exp  ^^Jdy 


(24) 


When  the  mass  flow  rate  is  zero  at  times  greater,  than  tqj  the  following  solutions 
for  f(T)  and  f'(T)  can  be  obtained. 
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f(T)  = f(T^)ex^  r.  ^ - T^)  1 (26) 

L o J 


Now  that  the  results  for  f(T)  and  f'(T)  are  available,  Eqs.  21  and 
22  for  the  overpressure  and  particle  velocity  cein  be  expressed  as  shown  below. 


if  T < 0 


Ap(t)  = < 


^ ^ exp  ( - -^  ) J [ffl»(y)N;(y) 
■"o  o 

+ ini(y)N^(y)]exp  ^ ^ ^ 


Ap(T^)e3q)  - ^ (t  - T^) 

L o 


if  0 < T < T 


if  T > T 


if  T < 0 


7^  exp  ^ r [ni^(y)N4(y)  + m;(y)N^(y) 
or  A \ /J 

o 

+ J m^(y)N^(y)]exi>  ^ ^ dy 


Au(t)  = ^ 


if  0 < T < T. 


Au(TQ)exp  - ^ (t  - T^)  if  T > 

L o J 


(28) 


(29) 


These  results  con5)lete  the  general  solution  for  the  wave  motion  in  the  pyramidal 
horn,  for  any  specified  time  distributions  of  the  mass  flow  rate  per  unit  area 
m^(T)  and  throat  area  N^(t) . In  most  cases  the  integrations  would  have  to  be 
done  numerically  to  obtain  final  results  for  the  overpressure  and  particle- 
velocity  signatures . 

Let  the  mass  flow  rate  per  unit  area  be  given  by  the  gasdynamic 
analysis  (Eqs.  15  and  l8) , and  also  let  the  throat-area  variation  be  parabolic 
(Eq.  12).  After  these  expressions  sire  substituted  into  Eqs.  28  and  29,  the 
integrals  can  be  evaluated  analytically.  The  integration  procedure  is  not 
unduly  difficult,  but  it  is  a rather  long  and  tedious  process.  The  final 
results  for  the  overpressure  signatiure  are  summarized  below,  in  a sinple  and 
convenient  form. 


19 


L 


£^(T)  = 0 


if  T < 0 


Ap(T)  = Zi[(l  - 2T  + 2p  - (1  + 2P)exp(-  T/P)} 
+ a^e[lOT^  - 20T^  + 9T^) 

+ a^€p(-40T^  + 60T^  - i8t) 

+ a^e^[-32I^  + 112T^  - 126t^  + 45T**) 

+ a^€p^(l20T^  - 120T  + l8(l  - exp(-  T/P))) 

+ a^e%[22kT^  - 672T^  + 6301^  - l80T^) 

+ a2€^{88T^°  - 4U0I^  + 8lOT®  - 648t^ 


189T°)] 


if  0 < T < 1 


Ap(T)  = Ap(l)exp[-(T  - l)/p] 


K 

o A 


if  T > 1 


T = — 

T 
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In  the  above  results  the  coefficients  a^^,  82  and  a^  are  not  new  but  have  been 
given  previously  in  Eq.  18,  and  ^ is  a characteristic  aa^jlitvide  of  the  over- 
pressure signature.  Note  that  both  parameters  P and  e are  typically  smAn 
compared  to  unity,  and  that  the  expression  for  the  overpressure  signature  is 
thus  accurate  to  third  order. 

The  characteristic  anplitude  ^ in  Eq.  30  can  be  expressed  in  an 
alternate  and  more  convenient  form.  Firstly,  A^Aq  can  be  replaced  by 
where  r^  is  an  equivalent  ra^us  corresponding  to  the  maximum  throat  area  A^. 
For  exasple,  if  A^  is  156  cm^  for  the  Travelling-Wave  Horn,  then  r„  equals 
1.0  m by  using  Eq.  1.  Secondly,  the  initial  mass  flow  rate  per  luiit  srea  iijj. 
or  m^(T  = 0)  can  be  eliminated  by  using  Eq.  15,  and  the  initial  Mach  nunfcer 
in  the  resxilting  expression  can  be  removed  by  using  Eq.  10.  Finally,  after 
some  rearrangement,  the  desired  expression  for  ^ is  given  below. 


/5P  = 


hJh.  r 

^ a_T.  i 


Hence,  this  res\ilt  is  expressed  in  terms  of  the  initial  reservoir  pressure  Po 
and  sound  ;^eed  Oq,  atmospheric  pressure  p,  flow  or  N-wave  duration  Tq  and 
equivalent  reuiius  r^.  ° 

The  expression  for  the  particle-velocity  signature  can  he  derived  in 
a manner  similar  to  that  used  to  obtain  the  ejQpression  for  the  overpressure 
signature.  However,  now  that  the  results  for  the  overpressure  are  known,  it  is 
much  sinpler  to  use  these  known  results  in  the  following  equation  to  obtain  the 
expression  for  the  particle-velocity  signature. 


T 

z^(T)  = r^(y)dy 

Q a or  d 


Because  the  integration  in  this  equation  can  be  done  rather  easily  to  obtain 
the  final  results  for  the  particle  velocity,  these  quite  lengthy  results  are 
not  presented.  Note  that  Eq.  32,  which  conveniently  expresses  the  particle 
velocity  in  terms  of  the  overpressure,  can  be  obtained  directly  from  Eqs.  21 
and  22. 

In  order  to  illustrate  some  interesting  and  inportant  features  of  the 
overpressure  and_p^ticle-velocity  signatures,  consider  the  following  sinplified 
case.  Let  e or  a^A^ro/V  be  zero,  which  corresponds  to  the  case  of  a reservoir 
having  an  infinite  volume.  Hence,  the  reservoir  conditions  such  as  the  pressure 
and  tenperature  are  constant  during  the  generation  of  a simulated  sonic  boom. 
Mso,  let  P or  ro/axo  be  zero.  This  corresponds  to  the  case  of  a point  source 
(Ref . 10)  or  the  reservoir  air  entering  the  horn  at  its  apex  (ro  = 0) . Then, 
the  res\ating  sinple  expressions  from  Eqs.  30,  31  and  32  for  the  overpressure 
and  particle-velocity  signatures  are  given  below. 


/^SpCI  - 2T) 


r ZPU 

I 0 


if  0 < T < 1 
if  T < 0 or  T > 1 


l-2T+~(T-T^)]ifO<T<l 

Au  = I p a L ^ 

0 ifT<0orT>l 


A nuntoer  of  inportant  observations  can  be  made  from  these  sinple  results. 
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The  wave  or  simulated  sc*iic  boom  that  propagates  from  the  cmai  i to  the 
large  end  of  the  pyramidal  horn  has  a predicted  N-shaped  signature  (Eq.  (33), 
which  is  illustrated  graphically  in  Fig.  jha.  Remember  that  the  throat-area 
variation  was  taken  to  be  parabolic  for  this  prediction.  The  two  in 

the  signature  have  an  instantaneous  or  aero  rise  time  and  an  identical  anplitude 
^ (Eq.  31) • This  an^litude  increases  with  increasing  reservoir-pressure  ratio 
Po/p,  and  it  is  inversely  proportional  to  both  the  flow  or  N-wave  duration  tq 
and  radius  Fq.  Hence,  the  overpressure  profile  is  invariant  with  distance  but 
the  overpressxire  diminishes  like  l/r. 

Althovigh  tte  particle-velocity  signature  also  exhibits  two  shocks 
of  equal  aii?)litude  (^/p  a)  and  instantaneous  rise  time,  it  is  not  necessarily 
N shaped  or  invariant  with  distance  (see  Eq.  34),  Velocity  signatures  for 

fixed  radii  are  shown  in  Fig.  34b.  All  of  these  signatures,  charac- 
terized by  different  values  of  the  parameter  aro/r,  have  been  normalized  and 
given  on  the  same  diagram  to  illustrate  the  change  in  wave  form  from  an  almost 
parabolic  shape  in  the  near- field  to  an  N-shaped  profile  in  the  far- field.  Often 
it  is  convenient  to  use  the  inverse  parameter  r/axQ  to  identify  a velocity  profile 
at  a relatively  large  radius.  Because  axo  equals  the  wave  length  >©,  the  inverse 
parameter  r/axo  ^/^o  is  a direct  measure  of  how  many  wavelengths  the  wave  is 
from  its  source. 

The  interior  test  section  of  the  UCIAS  Travelling-Wave  Horn  is  not 
located  in  the  far-field  but  rather  in  the  so-called  mid-field,  as  the  radial 
distance  to  this  section  (20  to  25  m)  is  only  a fraction  of  a wave  length.  For 
exaaple,  the  values  of  r/^  are  equal  to  0.45,  O.30  and  O.15  for  simulated  sonic 
booms  having  different  durations  of  100,  200  and  300  ms  respectively.  Conse- 
quently the  particle-velocity  profile  as  predicted  at  the  interior  test  section 
is  not  N sh^ed,  although  the  signatwe  is  not  too  different  than  the  N shape. 
Because  the  signature  has  a convex  curvature  between  the  two  shocks  the  positive 
velocity  portion  has  a greater  duration  than  the  negative  velocity  portion. 

These  resixlts  can  be  interpreted  from  the  velocity  profiles  given  in  Fig.  34b, 
by  noting  that  axo/r  equals  2,  3 and  6 for  the  respective  durations  of  100, 

200  and  300  ms.  These  particular  profiles  are  not  shown  for  clarity. 

More  interesting  and  iuportant  features  of  the  overpressure  signature 
can  be  illustrated  by  letting  3 or  rg/axo  zero  as  for  the  previous  case  but 
now  letting  e or  a^A^xo/V  be  nonzero.  Since  the  parameter  a„A„Xo/V  is  a func- 
tion of  the  reservoir  volume  (V),  N-wave  duration  (xo)  and  flap- valve  throat 
area  (A^) , their  influence  on  the  overpressure  signature  can  be  Jointly  evalua- 
ted. The  required  expression  for  the  overpressure  can  be  obtained  from  Eq.  30, 
and  it  is  given  below. 


/2p=0  if  T<0  or  T>1 
^ = ^[1  - 2T  + a^eClOT^  - 201^  + 9T^} 

+ a2€^(-32T^  + 112T^  - 126t^  + 45T^)  (35) 

+ a^e^CSSr^®  - 44ot^  + 81OT®  - 648t'^ 

+ 189T^)]  if  0 < T < 1 


22 


Ik 


Three  sets  of  graphical  illustrations  of  pressure  profiles  for  different 
typical  values  of  e or  ILa^to/V  of  0.05,  0.10  and  0.15  are  given  in  Fig.  35. 

In  each  set  of  results  the  dependence  of  the  profile  on  the  initial  reservoir 
pressure  ratio  po/p  is  illustrated.  It  can  be  seen  that  the  distortion  of  the 
averpressure_signature  from  the  ideal  N shape_becon»s  worse  with  increasing 
values  of  a^A^To/V  and  diminishing  values  of  Po/p.  These  results  illustrate 
that,  if  a long-duration  simulated  sonic  boom  having  a good  W- shaped  signature 
is  to  be  generated  with  a Travelling- Wave  Horn  utilizing  a large  value,  then 
the  reservoir  volpae  (v)  must  be  sufficiently  large  to  offset  large  values  of 
the  throat  area  (A^)  and  duration  (tq) • The  parameter  a^A^To/V  should  there- 
fore have  a value  which  is  less  than  about  0.1. 

The  results  shown  in  Fig.  35  have  been  rearranged  and  given  again 
in  Fig.  36.  The  three  sets  of  results  correj^ond  to  different  values  of  the 
reservoir  pressvire  ratio  po/p  of  1.007,  I.015  and  I.89.  In  each  set  the 
effects  of  different  values  of  a^A^To/V  of  0,  0.05,  0.10  and  0.15  on  the 
signature  are  shown. 

Further  interesting  and  inqportant  features  of  the  oveipressure 
signat\jre  can  be  illustrated  by  letting  P or  rg/aTo  nonzero  and  e or 
a^A^To/V  be  zero  (infinite  reservoir  volume).  For  this  case  the  following 
expression  for  the  overpressure  can  be  obtained  frcai  Eq.  30. 


Ap  = 0 if  T < 0 

AP  = - 2T  + 2p  - (1  + 2P)e:5)(-T/p)]  if  0 < T < 1 (36) 

^ = ^[-1  + 2p  - (1  - 2p)exp(-l/p)  exp  -(T  - l)/P]  if  T > 1 

Based  on  these  results  the  effects  of  a nonzero  value  of  p or  rg/aro  on  the 
N- shaped  overpressure  profile  can  be  assessed,  and  these  effects  are  conveniently 
ill^trated  graphically  in  Fig.  37*  The  essential  result  of  a nonzero  value  of 
rg/axg  is  that  the  front  and  rear  shocks  of  the  profile  both  have  a finite  rise 
time.  It  can  be  seen  that  the  two  rise  times  become  longer  relative  to  the 
N-wave  duration  as  the  value  of  rg/axg  increases.  Actually  the  rise  time  is 
nearly  constant,  but  it  aTOears  more  conspicuously  on  shorter  duration  N-waves 
for  which  the  vilue  of  ro/axg  is  larger.  Note  that  an  instantaneous  rise  time 
would  occur  only  when  the  value  of  rg/axg  or  rg  is  zero.  (The  radius  rg  is 
the  location  at  which  the  mass-flow  valve  is  Joined  to  the  pyramidal  horn.) 

The  peak  overpressure  of  the  front  shock  of  an  N-wave  having  a finite 
rise  time  is  always  less  than  that  having  an  instantaneous  rise  time,  and  this 
difference  becomes  more  marked  for  a shorter  duration  N-wave.  This  featvure 
is  also  illustrated  by  the  three  overpressure  signatures  shown  in  Fig.  37. 
Quantitatively,  the  maxi mum  overpressure  of  the  front  shock  and  the  corre- 
sponding rise  time  At®  from  zero  to  the  maximxan  overpressure  can  be  derived 
easily  from  Eq.  36,  and  these  two  results  are  given  below. 
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These  results  demonstrate  that  the  rise  time  Atjq  increases  and  the  peak  over- 
pressure ^>111  diminishes  as  the  N-wave  duration  tq  increases. 


Previously,  in  Chapter  3»  ‘tke  rise  time  of  a shock  was  defined  to  be 
equal  to  1.25  times  the  time  for  the  overpressvire  to  rise  from  10^  to  ^Oio  of 
its  peak  value.  By  using  Eq.  36,  the  respective  times  ‘'^90 

overpressure  of  the  front  shock  to  rise  frctu  zero  to  10%  and  90%  oT  the  peak 
value  are  given  by  the  following  approximate  expressions. 


90 


■^10- 
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1 - 2P  in 


1 + 2P 


0.1  + 2P  + l.l6^  in 


1 + 2P 

f \ + 2&  \ - 


These  resialts  are  accvirate  for  values  of  p or  Tq/^tq  ranging  from  0 to  0.2, 
which  are  more  than  sufficient  for  the  present  work.  The  rise  time  Axp  or 
1.25(t^O  “ Tio)  therefore  given  by  the  following  approximate  expression. 


At„  « 1 
r 


1 + 2p 

0.1  + 2p  + I.16P  in 


1 + 2P  ^ ^ 


-0.1  + 0.2P  in 


1 + 2P 


Although  the  rise  time  Arp  diminishes  monotonically  with  increasing  val\xes  of 
P,  Atp  is  essentially  constant  at  2.75  ro/a  for  small  P values  (<  0.005)  • 
Furthermore,  for  P values  ranging  from  0.01  to  0.05,  which  are  typical  for  the 
l/EIAS  Travelling-Wave  Horn,  Arp  is  roughly  equal  to  2ro/a.  Results  such  as 
these  can  be  determined  from  Eq.  39>  or  from  Fig.  38  where  Arpi/rp  is  plotted 
versus  P or  ro/arQ.  Note  that  for  rg  and  a values  of  I.58  m and  430  m/s 
respectively,  the  predicted  rise  time  is  roughly  9 uis. 

As  a final  exanple  of  the  analytical  results  for  the  overpressure 
signatxjre  of  a simulated  sonic  boom,  three  typical  signatures  for  the  ITPIAS 
Travelling-Wave  Horn  are  shown  in  Fig.  39 • It  can  be  seen  that  the  three 
N-shaped  profiles  (solid  lines)  have  a nominal  peak  overpressure  of  100  N/n^ 
and  different  durations  of'  100,  200  and  300  ms.  They  exhibit  finite  rise 
times  (P  ji^  0)  and  distortions  from  the  N shape  due  to  a finite  reservoir 
volume  (e  ^ 0) . These  slgnatures_were  calcvilated  by  using  Eq.  30  and  typical 
values  of  p or  Tq/b.tq  and  e or  a^A^Tg/V.  In  eacn  (.ase  the  corresponding  ideal 
N-wave  (dotted  line)  featuring  an  instantaneous  rise  time  (P  = 0)  and  no 
distortion  (c  = 0)  is  shown  for  completeness  and  comparison  pxirposes. 

4.3  Comparison  of  Predicted  and  Measured  Results 

All  of  the  predicted  overpressure  signatures  which  are  compeired  to 
experimental  results  in  this  section  have  been  calciilated  by  using  Eqs.  30  and 


31.  Iii5)licit  in  the  analysis  and  thus  these  equations  is  the  assuiuption  that 
the  throat-area  variation  with  time  of  the  flap  valve  is  parabolic,  as  given 
by  Eq.  12.  The  maximum  throat  area  for  the  IfEIAS  Travelling- Wave  Horn  and 
the  corresponding  equivalent  radius  r^  used  in  the  calculations  were  156  cm2 
and  1.0  m respectively,  and  the  radial  location  ro  at  which  the  valve  is 
Joined  to  the  pyramidal  horn  was  taken  to  be  1.58  m.  Italess  noted  to  the 
contrary  the  reservoir  voliane  V was  10.2  m3.  Furthermore,  the  overpressxare 
signatures  were  predicted  for  a radial  location  r of  21.3  in,  corresponding  to 
the  interior  test  section  of  the  horn  where  the  overpressure  measurements 
were  made. 

Three  measured  overpressure  signatures  having  the  same  nominal 
duration  of  100  ms  have  been  traced  from  the  oscillograms  in  Fig.  2k  and 
reproduced  in  Fig.  4o.  Corresponding  predicted  signatures  are  also  shown  in 
Fig.  40,  in  order  to  facilitate  their  comparison  to  measiired  profiles.  Two 
similar  sets  of  measured  signatures  traced  from  oscillograms  of  Figs.  25  and 
26  for  nominal  dxirations  of  200  and  300  ms  are  reproduced  in  Figs.  4l  and  42 
respectively,  where  the  predicted  signatures  are  also  given.  It  can  readily 
be  seen  in  all  three  sets  of  results  that  the  predicted  and  measured  signatures 
are  in  good  agreement.  In  general  the  overpressure  signatures  measured  in  the 
UCIAS  Travelling-Wave  Horn  can  be  predicted  successfully  by  the  analysis 
developed  previously  in  this  chapter. 

Two  additional  measured  oveipressure  signatures  having  a nominal 
duration  of  200  ms  are  shown  in  Fig.  43.  For  this  case  the  reservoir  voliane 
of  the  facility  was  reduced  by  one- third  from  10.2  to  3*4  m^.  Included  in 
the  figiire  are  the  corresponding  predicted  profiles,  which  are  again  in  good 
agreement  with  the  measxired  signatures.  Note  that  both  the  measured  and  predicted 
results  differ  fairly  markedly  from  the  desired  N shape  for  good  simialation  of  a 
sonic  boom,  owing  to  an  inadequate  reservoir  volume.  Such  results  illustrate 
the  importance  of  having  a reservoir  with  sufficient  voliane. 


5.  DISCUSSIONS  AND  CONCLUSIONS 

The  UTIAS  Travelling-Wave  Horn  has  proven  to  be  a practical  facility 
for  the  simialation  of  a full-scale  sonic  boom.  The  air  compressor  and  reservoir 
system,  flap  valve  and  electric-drive  system  (including  the  flywheel,  clutch, 
cam,  three-bar  linkage,  two  flaps  and  central  fairing),  pyramidal  horn  and 
porous-piston  type  of  reflection  eliminator  all  essentially  function  correctly 
as  designed.  The  simulator  therefore  has  the  desired  capability  of  producing  a 
good  N-shaped  overpressiire  signature  or  simulated  sonic  boom  in  the  interior 
test  section  of  the  pyramidal  horn.  The  N-wave  amplitude  and  duration  can  be 
controlled  individually  to  be  less  than,  equal  to,  or  more  than  that  of  an  actual 
sonic  boom.  Owing  to  certain  design  features  of  the  compressor,  reservoir, 
flap  valve  and  control  system,  the  facility  has  the  desirable  capability  of 
producing  respectable  and  repetitive  simulated  sonic  booms  (l  to  10  booms  per 
minute),  in  order  to  facilitate  structural  response,  fatigue  and  damage  studies. 
Note  that  the  primary  use  of  the  ITEIAS  Travelling- Wave  Horn  is  to  study  the 
effects  of  simulated  sonic  booms  on  structiures  and  their  subcomponents.  Hence, 
the  generation  of  repeatable  and  repetitive  N-waves  at  a high  repetition  rate  is 
a necessary  reqviirement  of  the  facility. 


One  disadvantage  of  the  UTIAS  Travelling-Wave  Horn  utilizing  the  present 
flap  valve  is  that  the  rise  times  of  the  front  and  rear  shocks  of  the  simulated 


scaiic  boom  axe  longer  at  6 to  10  ms  than  those  of  an  actual  sonic  boom  (typically 
about  1 ms) . Also,  the  rise  times  are  not  easily  controlled  but  are  essentially 
fixed  by  the  valve  design  and  N-wave  duration.  Since  the  rise  times  normally 
have  very  little  effect  on  structural  response,  the  major  vise  of  the  facility 
for  structural  response,  fatigue  and  damage  studies  is  not  inpaired  by  the  present 
long  rise  times  of  6 to  10  ms,  or  by  not  being  able  to  control  the  rise  times. 
However,  human  and  animal  response  tests  requiring  a simulated  sonic  boom  having 
a short  rise  time  which  can  also  be  varied  from  0.5  to  10  ms,  for  exasple,  cannot 
be  done  with  the  present  facility.  However,  with  additional  developmental  costs 
the  rise  times  could  be  inproved  if  the  need  should  arise.  (Note  that  such 
response  tests  can  be  done  with  the  cooplementary  IHIIAS  Loudspeaker-Driven  Booth.) 

It  is  worth  mentioning  that  the  original  flap  valve  utilired  a specially 
designed  central  faiiring  (Fig.  27),  which  resulted  in  a short  rise  time  as  low  as 
1 ms  for  the  front  shock  of  the  simulated  sonic  boom.  Furthermore,  this  central 
fairing  haul  special  design  features  which  permitted  the  rise  time  to  be  controlled 
from  1 to  15  ms.  However,  because  this  central  fairing  suffered  frequent  mechani- 
cal failure,  a newer  present  centraJ.  fairing  (Fig.  l8)  having  good  mechanical 
reliability  was  designed  and  used  in  the  flap  valve.  This  new  centrail  fairing 
resulted  in  the  undesirably  longer  rise  times  of  6 to  10  ms.  Hence,  short  rise- 
time capability  wan  sacrificed  to  obtain  a nonfalling  mechanical  operation.  It 
is  believed,  however,  that  such  a sacrifice  is  not  necessary.  A future  modifi- 
cation to  tile  central  fairing  or  a whole  new  design  should  result  in  both  good 
mechanical  reliability  and  a short  rise  time  control  of  the  rise  time  from  1 to 
15  ms,  for  exaaple,  should  also  be  possible. 

The  Travelling-Wave  Horn  utilizing  a large  air  reservoir  and  large- 
^ea  flap  valve  produces  a simulated  sonic  boom  with  very  little  svperposed 
jet  noise.  Hence,  the  main  problem  of  intense  jet  noise  svperposed  on  the 
N-wave,  when  the  facility  used  the  previous  small-area  plug  valve  (Ref.  10), 
has  been  virtually  solved.  By  virtue  of  the  large- area  flap  valve,  the  mass- 
flow  rate  of  air  required  to  produce  a simulated  full-scale  sonic  boom  can  be 
delivered  at  a very  low  flow  velocity.  The  low  flow  velocity  at  the  flap  valve 
is  the  main  reason  that  insignificant  jet  noise  is  produced.  It  should  also  be 
noted  here  that  the  use  of  a large  valve  results  in  the  use  of  a low  reservoir 
pressure.  However,  this  advantage  is  offset  because  the  reservoir  has  to  have 
a large  air- storage  capacity,  such  that  the  reservoir  pressvire  does  not  change 
significantly  during  the  generation  of  a simulated  full-scale  sonic  boom. 

The  combined  gasdynamic  and  acoustic  analysis  successfully  describes 
the  wave  motion  in  the  pyramidal  horn  during  the  generation  of  a simvilated  sonic 
boom.  On  the  one  hand  the  gasdynamic  analysis  has  the  capability  of  predicting 
the  pressure,  terperature  and  other  state  properties  of  the  ejjpanding  reservoir 
air  and  ^so  the  mass  flow  rate  and  flow  Mach  nvimber  of  the  air  passing  through 
the  opening  and  closing  flap  valve.  On  the  other  hand  the  acoustic  analysis 
is  able  to  predict  the  overpressure  and  particle- velocity  signatures  of  the 
simulated  saaic  boom  in  the  pyramidal  hern,  including  their  main  characteristics 
such  as  the  peak  overpressure,  duration  and  rise  time.  One  inportant  result 
of  the  analysis  is  a suitable  criterion  on  how  large  the  reservoir  voliane  must 
be,  in  order  that  the  overpressure  signature  is  not  significantly  distorted 
from  the  N shape.  For  normal  operation  of  the  UTIAS  Travelling- Wave  Horn  the 
nondimensional  parameter  a^A^To/V  should  be  less  than  0.05,  0.10  and  0.15  to 
avoid  significant  distortion  in  an  N-wave  having  a duration  of  100,  200  and 
300  ms,  respectively.  Note  that  the  respective  symbols  a.^,  A^,  Tq  and  V denote 
sound  speed  at  the  valve  throat,  maximum  valve- throat  area,  N-wave  duration  and 
reservoir  volume. 
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FIG.  3.  ELEVATION  AND  PLAN  VIEWS  SHOWING  MOST  OF  THE  RELEVANT 

ELEMENTS  OF  THE  TRAVELLING-WAVE  HORN  USING  A PLUG  VALVE. 
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SIMPLIFIED  ASSEMBLY  DRAWING  OF  THE  PLUG  VALVE  EXHIBITING 
THE  OPENING  (a)  AND  CLOSING  (b)  MOVEMENTS. 
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SIMULATED  SONIC  BOOMS  GENERATED  WITH  AND  WITHOUT  THE 
JET-NOISE  ABSORBER  (plug  valve,  radial  distance  of  15.2  m 
reservoir  pressures  of  1.3,  2.5,  5.1,  9.2  and  13.9  atm). 
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PLAN  VIEW 


FIG.  9. 


PLAN  VIEW  ILLUSTRATING  THE  MOST  RELEVANT  ELEMENTS  OF  THE 
TRAVELLING-WAVE  HORN  USING  A FLAP  VALVE. 


FIG.  10.  OUTSIDE  VIEW  OF  THE  TRAVELLING- WAVE  HORN  SHOWING  THE 

CONCRETE  HORN,  TEST  AND  CONTROL  ROOMS,  COMPRESSOR  SHED, 
AND  TWO  TANKS  OF  THE  AIR  RESERVOIR  SYSTEM. 


FIG.  11.  VIEW  O'"  THE  SMALL  END  OF  THE  HORN  SHOWING  ALSO  THE  LARGE  TANK 
OF  THE  AIR  RESERVOIR  SYSTEM,  FLAP  VALVE  AND  ELECTRIC-MOTOR 
DRIVE  SYSTEM,  CONTROL  CONSOLE,  AND  ANCILLIARY  ELECTRONICS. 


FIG.  12.  VIEW  OF  THE  LARGE  END  OF  THE  HORN  SHOWING  ALSO  THE 
REFLECTION  ELIMINATOR,  ACCESS  DOOR  TO  THE  HORN 
INTERIOR,  AND  CUTOUT  FOR  STRUCTURAL  PANEL  TESTS. 
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RESERVOIR 


FIG.  15.  FLAP  VALVE  OPERATION:  a)  VALVE,  b)  THROAT-AREA  AND 
MASS-FLOW-RATE  PROFILES,  c)  OVERPRESSURE  SIGNATURE. 


TRIGGER  MECHANISM 


FIG.  22.  MEASURED  ANGULAR-DISPLACEMENT  HISTORIES  OF  A VALVE 
FLAP  FOR  THREE  DIFFERENT  INDICATED  DURATIONS  AND 
RESERVOIR  PRESSURES. 
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FIG.  24.  SIMULATED  SONIC  BOOMS  HAVING  A DURATION  OF  100  ms 
(measured  at  a radial  distance  of  21.4  m). 


FIG.  25.  SIMULATED  SONIC  BOOMS  HAVING  A DURATION  OF  200  ms 
(measured  at  a radial  distance  of  21.4  m). 
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FIG.  26.  SIMULATED  SONIC  BOOMS  HAVING  A DURATION  OF  300  ms 
(measured  at  a radial  distance  of  21.4  m). 


FAIRING  PLATE 


FIG.  27. 


MOTION  OF  THE  CENTRAL-FAIRING  PLATES  AND  VALVE  FLAPS 
FOR  A CENTRAL-FAIRING  DESIGN  WHICH  GAVE  A FAST  RISE 
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FIG.  40.  MEASURED  AND  PREDICTED  OVERPRESSURE  SIGNATURES  HAVING 
A NOMINAL  DURATION  OF  100  ms  (A*=  156  cm  , r*=  1.0  m, 
= 1.58  m,  r = 21.3  m,  V ^ 10.2  m*,  a = 340  m/s). 


